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Abstract - We study a process termed agglomerative percolation (AP) in two dimensions. In- 
stead of adding sites or bonds at random, in AP randomly chosen clusters are linked to all their 
neighbors. As a result the growth process involves a diverging length scale near a critical point. 
Picking target clusters with probability proportional to their mass leads to a runaway compact 
cluster. Choosing all clusters equally leads to a continuous transition in a new universality class 
for the square lattice, while the transition on the triangular lattice has the same critical exponents 
as ordinary percolation - violating blatantly the basic notion of universality. 
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Percolation is a pervasive concept in statistical physics 
and an important branch of mathematics [1]. It typifies 
the emergence of long range connectivity in many systems 
such as the flow of liquids through porous media [5] , trans- 
port in disordered media [3], spread of disease in popula- 
tions [1] , resilience of networks to attack [5] , formation of 
gels [6] and even of social groups [7]. It also underlies a 
number of other critical phenomena - like the Ising or- 
der/disorder transition, which is a percolation transition 
on the set of spins with given sign [5] . 

■ The phase transition in ordinary percolation (OP), 
where bonds or sites are added at random, represents a 
broad universality class. Recently Achlioptas et al. [3] 
rnade a simple modification by, at each step, selecting 
among two possibilities the link that leads to the slowest 
growth of large clusters. This global choice introduces a 
large length scale - the system size L - which can alter uni- 
versality. Indeed, they concluded that an unusual, discon- 
tinuous transition (called "explosive percolation" ) emerges 
where a macroscopic cluster appears suddenly while at the 
same time scaling in other quantities is observed |10H12) . 
Various modifications of the rule have been made P^IT^ 
- all finding evidence of a discontinuous transition. Al- 
though the claim for discontinuity in [5] was refuted later 
in |15j , explosive percolation in that case does represent a 
new universality class. 

Here we discuss a percolation process that also contains 
a potentially large length scale in its definition, in this case 



the correlation length ^. Our process has direct applica- 
tion to the study of complex networks. Instead of adding 
bonds randomly, we pick a random cluster and add bonds 
to its entire surface in order to link it to all adjacent clus- 
ters. Starting with the state where all clusters have size 
one, at each update t ^ t + 1 the process repeats until 
the entire lattice (or graph) is reduced to a single cluster. 
We call this "agglomerative percolation" (AP), in anal- 
ogy with cluster growth by aggregation [TB]. Thus if by 
chance a cluster of length scale i is picked, links are added 
simultaneously at distances 0{l) apart. 

AP can be analyzed on any graph. It corresponds 
to random sequential renormalization jl7| of a network, 
where a single cluster is identified as a 'super'-node that 
is a local coarse-graining of the graph. In this perspec- 
tive, scaling laws seen in renormalization studies of small- 
world networks [TBHH] are a consequence of an AP phase 
transition and do not indicate fractality of the underly- 
ing graph ^Jl^^. Previously AP was studied on critical 
trees [17j and in one dimension |23j . Scaling laws were 
found both analytically and numerically - but no phase 
transition occurs since both graphs have a topological di- 
mension of one. 

In order to establish the phase transition in AP, and its 
relationship to OP, we analyze it in two dimensions (which 
is clearly not a fractal graph), where many exact results 
for OP are known. We consider both square and triangu- 
lar lattices. Clusters can be chosen with equal probability. 
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Fig. 1: (Color online) The six largest clusters in mass- weighted 
AP (model (b)) on a square lattice with L — 256 (top), and 
L — 8192 (bottom). The average cluster mass is two, or n = 
N/LF' = 0.5. The red cluster is close to wrapping. All large 
clusters are compact. 

or we can make biased choices according to the mass, ra- 
dius, etc. of the clusters. Here both (a) uniform probabil- 
ities and (b) probabilities proportional to the cluster mass 
are studied. Model (b) coincides with choosing sites uni- 
formly, and growing the whole cluster in which they lie. It 
shows runaway behavior resembling a first order transition 
with compact clusters - as in Ref. [23] (see Fig. 1). Model 
(a) is more subtle. Although clusters appear fractal (see 
Fig. 2) and the overall character resembles OP, fundamen- 
tal differences arise. Most conspicuous is an unexpected 
difference between the two lattice types: While model (a) 
on the square lattice is definitely not in the OP univer- 
sality class (e.g. the average cluster size diverges at the 
transition), the triangular lattice shares the same critical 
exponents as OP. We believe this violation of universality 
must be related to the long range nature of the growth 
process for large clusters. 

We use L X L lattices, with 2^ < L < 2^^. Boundary 
conditions are helical: periodic in the y direction while the 
right neighbor of site (x = L, y) is (x' = 1, = ?/(mod L)-|- 
1). Diagonal bonds are added to obtain triangular lattices. 
We use an algorithm based on that in [53] , augmented by a 
depth first search on the target cluster, in order to find all 
its neighboring clusters. The natural control parameter 



in OP is p, the fraction of existing bonds or sites. In 
AP, however, the number of links is not uniquely defined, 
although it is a version of (correlated) bond percolation. 
One might join two clusters via a single link, but one might 
also put multiple links between them. Therefore, in AP 
it is more natural to use the average cluster number per 
site, 

n = N / ^ {m)'^ , (1) 

where (to) is the average cluster mass. For OP, n{p) is not 
analytic at p = Pc, but is monotonic with two continuous 
derivatives. Thus one can use n as the control parameter 
in OP and reproduce all known scaling laws. (We checked 
this explicitly; see also [lU])- Instead of n, one might also 
use the number of agglomeration events, t, as the control 
parameter. In agreement with |17[|26| . we found that t is 
more "noisy" than n and leads to slightly less clear results. 
We analyze the distribution Pn{m) of cluster masses to in 
configurations with cluster density n, and the probability 
Pwrap('T-), that a cluster wraps the torus in the j/— direction, 
for each L. 

In model (b), the growth rate for a cluster of mass m 
accelerates steeply with to, leading to a runaway effect. 
A cluster's chance to be selected is mjl? . Once chosen 
it grows all along its perimeter. Since most of its neigh- 
bors are small, it grows into a compact shape. For any 
< a < 2 and for i — oo, we conjecture that the largest 
cluster reaches mass to ~ L" at a time when (to) — ^ 1. 
This "incipient" cluster continues to separate in mass from 
the others. It wraps the torus when mjl?' = 0(1). Thus 
an infinite incipient cluster appears at density ric —> 1 in 
the limit L — > oo, while wrapping occurs much later, at 
< ric < 1. Fig. 1 shows the six largest clusters in a 
typical run on a square lattice for both small and large 
L at n = 0.5. These snapshot were taken at a time long 
past the appearance of the incipient cluster and long be- 
fore it wraps. One sees that the giant cluster becomes 
more dominant over all other clusters as L increases. Al- 
though convergence of ric 1 as L —> oo is slow, it is in 
perfect agreement with numerical simulation results (data 
not shown). The same scenario holds for the triangular 
lattice. 

In OP, cluster perimeters are for large clusters propor- 
tional to their mass. Thus, if a new bond is added at 
each time step, the average growth rate of a cluster is 
roughly dm/dt oc to. In AP model (a) - where clusters 
are picked with uniform probability - those chosen grow 
by an amount proportional to their perimeter, so again 
(roughly) dm/dt oc to. This leads neither to a runaway 
of large clusters as in model (b) nor to the retardation 
of their growth as in the Achlioptas process. Therefore, 
model (a) and OP cannot be distinguished by such a crude 
argument and their relationship could conceivably depend 
on microscopic details such as the type of lattice. 

For model (a), we first consider square lattices. By eye 
individual configurations look like OP. However, wrapping 
thresholds depend strongly on L. Figure 2 displays the 
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Fig. 2: (Color online) The six largest clusters in typical runs 
of AP model (a) on square lattices for n = N/L^ = 0.1 — 
l/(m). Top panel: for L = 256, the red cluster has already 
wrapped; Bottom panel: for L — 16384, all clusters are far 
from wrapping. In both cases, clusters appear fractal. 

six largest clusters in a typical run when n = 0.1 for both 
L = 256 (top) and for L = 16384 (bottom). While the 
largest cluster clearly wraps the small lattice, it is far from 
this point on the large one. Figure 3 shows the density 
nc,wTa.p{L) at which half the runs contain a wrapping clus- 
ter. The data fall roughly on a straight line on a log-log 
plot. If deviations from a straight line were typical finite 
size corrections, this would mean that the average cluster 
size at the wrapping threshold diverges as a power of L 
when L — > oo. However, for reasons explained below, we 
believe that nc,map{L) — > logarithmically as i — >• oo (for 
explicit fits, see the supplementary material [17]). This 
implies that the correlation length exhibits an essential 
singularity as n ^> ric = 0. 

Mass distributions Pnini) for n ~ ?^c.wrap(i) are dis- 
played in Fig. 4 using a data collapse method which com- 
pares m7Pn(wL) to m/L'^ , with r = 1.926 and D = 1.96. 
A perfect collapse corresponds to a finite size scaling (FSS) 
ansatz 

P„im)^m-^f{^b{n,L),m/L'') , (2) 

which generalizes the standard FSS ansatz [5] where 
il){n, L) = (n — nc)L^^^ ■ Except for peak heights the data 
collapse is excellent. The apparent values for r and D 
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Fig. 3: (Color online) Critical densities nc{L) vs. L for AP 
model (a) on square lattices. The upper curve is the average 
wrapping threshold. The lower one is the density at which 
the power law range in the cluster mass distribution extends 
furthest. Error bars are smaller than the symbol size. 

deviate significantly from their values in OP (r = 2.055 
and D = 1.89 in two dimensions). At n = ?^c,wrap(i), 
Pn{m) is not horizontal over a wide range of masses. 
Hence n = nc,wiap(i) is not equal to ncriL). The lat- 
ter is the density at which a power law in Pn{m) extends 
over the broadest range. Curves for P„-„^ ^'"^ 
also shown in Fig. 4 and exhibit data collapse with a re- 
markably wide power law range. The power law regime 
describes the relatively few remaining clusters in config- 
urations dominated by one wrapping cluster. The values 
ric^riL) are also shown in Fig. 3. and decrease similarly to 
"■c,wrap(i) as L increases. 

If T < 2 as suggested by Fig. 4, then the average cluster 
size at criticality diverges as i — ^ oo, in agreement with 
Fig. 3 but in stark contrast to OP. Accepting this, the 
two possible (scaling) scenarios are r = 2 or r < 2. If 
T < 2, nc{L) vanishes as ~ with some 5 > 0, and 
apparent values for t and D should not vary much with 
L. Neither of these statements is correct. Figure 3 shows 
definite curvature, and the best fit values for r and D both 
increase slightly but significantly with L, see Fig. [S] and 
P7]). While these small shifts arc not visible on the scales 
shown in Fig. 4, they do not diminish as L increases. One 
would not expect to see large corrections to scaling that 
could explain Fig. [5] if r > 2 since in that case the average 
cluster size is finite. 

Since the numerical value of D is determined from the 
positions of the peaks in Fig. 4, D is actually the fractal 
dimension of the largest cluster. The contribution to (to) 
from this cluster is Sma^/N — s maxim) / L"^ . If one assumes 
Eq. (HI) and r < 2, then one can show that the largest 
cluster makes a non-vanishing contribution to (m) as L — >■ 
oo. This can only happen if Smax ^ L^, showing that the 
largest cluster has D = 2 if r < 2. Furthermore, one would 
expect convergence of the apparent D to follow a power 
law in that case. But the slow convergence of D from 
below indicates again that the behavior is dominated by 
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Fig. 4: (Color online) Data collapse for AP model (a) on square 
lattices: m'^ P„(m) vs. m/L° with r = 1.926 and D = 1.96. 
The curves with the smaller peaks are for n = nc,wrap(i), while 
the others are for n = nc,T{L). The first are at the wrapping 
threshold, while the second are when Pn(m) has the broad- 
est power law range. The straight tilted line corresponds to 
Pn{m) oc m^^. System sizes are L — 256, 512, . . . 16384. 

logarithms. 

Wc conclude thus that the true asymptotic values are 
T = D = 2 and v = oo. In addition, we measured the 
exponent a In agreement with the scaling relation 

(7 = D/iy, we found a ~ with rather slow convergence. 
Thus, all scaling relations are (trivially) satisfied (notice 
that the other exponents give no constraint, in the present 
case, on the order parameter exponent (3). 

For the triangular lattice, clusters look like those in OP. 
Both ric^riL) and nc,wrap(-^) converge rapidly to the same 
(finite) critical value ric = 0.1561 ± 0.0002. Indeed, the 
best estimates of t{L), obtained by fitting power laws to 
Pn=n^ also converge rapidly to r = 2.057±0.002, 

in perfect agreement with OP (see Fig. [5]). Also D, ob- 
tained from a data collapse as in Fig. 4, and the exponent 
a are both within error equal to their values in OP, al- 
though these error bars are larger than for r, see [27] . 

Small apparent inconsistencies with OP arise when we 
try to estimate v using the scaling hypotheses ric^riL) — 
Tic L"^/" or Pmap{n) = cj)[{n — nc)L^^'']. The first rela- 
tion gives V = 1.10 ± 0.07, significantly smaller than the 
value V = 4/3 for OP. The second one gives ly = 1.47±0.05 
for n > ric- Taken at face value, these estimates would ex- 
clude universality with OP. But we believe that they are 
artifacts of large finite size corrections. Figure |6] shows 
an attempted data collapse for dp^ra.p{n)/dn. While the 
collapse is acceptable for n > ric, huge tails develop for 
n ^ ric as L increases. For small n these tails decay 
roughly as [{ric ~ n)L°-^]~^ with fi « 1.5. Wc checked 
explicitly that the tails result from events where wrapping 
happened when a large target cluster was hit. In such 
cases, TV can make huge jumps, so that the largest n at 
which the cluster wraps is far below the actual threshold. 
Thus we conclude that AP on the triangular lattice is in 
the OP universality class, with the caveat that the def- 
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Fig. 5: (Color online) Plots of refi(L), the effective Fisher ex- 
ponents estimated from the longest stretches in P„(m) that 
are compatible with pure power laws. As these estimates are 
somewhat subjective, the error bars are subjective as well. But 
their order of magnitude is consistent with the smoothness of 
the data with varying L. Clearly there is non-trivial L depen- 
dence both for the triangular and for the square lattice, with 
Tefi(i) increasing with L in both cases. While Tcff is compat- 
ible with the value r = 187/91 = 2.0549 ... in the case of the 
triangular lattice, it is much smaller for the square lattice. 

inition of n is affected by occasional large jumps which 
do not modify the main critical exponents but which do 
modify the tails of scaling functions as in Fig. [HI 

For an appreciation how different the behaviors are on 
the square and triangular lattices, we show in Fig. [7] effec- 
tive critical cluster densities. We plot them against L~'^/*, 
since this should give for OP straight lines, according to 
the FSS ansatz. While the data for the triangular lat- 
tice indeed follow roughly straight lines and give a finite 
non-zero value of ric, the same is definitely not true for 
the square lattice: Those curves strongly bend down for 
L ^ oo, suggesting that ric = l/i^ = 0. Whether the lat- 
ter is correct or not, this figure should leave no doubt that 
the square lattice model is not in the ordinary percolation 
universality class. 

In summary, we have studied agglomerative percolation 
(AP) in two dimensions. This class of models is equivalent 
to random sequential renormalization schemes |17j first 
introduced to scrutinize the supposed fractality [THH12] of 
real world - in particular, small- world ~ networks. Regular 
lattices were chosen for two reasons: (1) They are not 
fractal; (2) Detailed comparison can be made with exact 
results for ordinary percolation (OP). Our results display 
some of the rich behavior possible in this general class 
of models and indicate that at least some of the scaling 
behavior found in [T8H22] is due to AP rather than any 
supposed fractality of the underlying graph. If clusters are 
chosen with a bias for larger mass (model (b)), a runaway 
effect separates the largest cluster from the others and 
the behavior is completely different from OP. If clusters 
are chosen with equal probability (model (a)), then only 
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Fig. 6: (Color online) Data collapse plot for the probability 
density dp^ra.p{n) / dn against (n — nc)L^^" with ric — 0.1561 
and u = 1.47 for model (a) on triangular lattices. The scale on 
the y-axis is adjusted such that all curves peak at y = 1. 

a detailed numerical scaling analysis shows that AP is not 
in the OP universality class on the square lattice. On the 
other hand, AP on the triangular lattice shares critical 
exponents with OP. 

AP may have applications beyond network renormaliza- 
tion. Growing clusters appear in many different physical 
situations. It could happen that further growth is trig- 
gered by some excitation where the entire cluster suddenly 
invades neighboring clusters at its boundary. Agglomera- 
tive percolation could also describe the growth of countries 
or urban areas. Countries often grow by overrunning and 
incorporating neighbors during aggressive periods, when 
they attack and incorporate simultaneously several of their 
neighbors. 

Acknowledgement: We thank Bob ZifF for correspon- 
dence that helped us sharpen some of our arguments. 
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In this supplementary material, we present several plots. They illustrate claims for which in the main paper either 
the actual data were not shown or they were plotted in different ways. The present plots are to convince the reader 
that: 

• Model (b) (where clusters are selected for agglomeration according to their mass) develops runaway clusters 
very early in the process, leading to a "phase transition" at (m) = 1 in the infinite system limit N oo (here, 
(to) is the average mass per cluster) (Fig. SI); 

• The same is not true for model (a), but model (a) shows qualitatively different behaviors on the square and 
triangular lattices (Figs. S2 to S5); 

• In particular, while model (a) on the triangular lattice is in the ordinary percolation universality class, it is 
definitely not in that universality class on the square lattice; 

• For the square lattice, superficial analyses would give a Fisher exponent r and a fractal dimension D that are 
both smaller than 2, but a more careful analysis, combined with analytic arguments, leads to the conclusion 
that T ~ D = 2 exactly. 
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FIG. 1. (Color online) Mass distributions for model (b) with L = 4096. Each curve corresponds to a fixed value of n (number 
of clusters per site), corresponding to an average cluster mass (m) — 1/n. The uppermost (red) curve shows a clear narrow 
peak, corresponding to a single cluster that involves a substantial (and thus not strongly fluctuating) fraction of all sites. The 
middle (green) curve shows a broad bump at its upper end, corresponding still to a runaway cluster in most realizations, albeit 
with a much wider mass range centered typically at a few percent of the total mass. The lowest (blue) curve corresponds to 
(m) « 1.14, i.e. by far most clusters consist of a single site. In spite of this, it shows a long tail implying that some clusters 
have reached size 10'' and beyond. In addition to Fig.l in the main paper, this shows most directly the existence of runaway 
clusters at values of (m) that converge to 1 as A'' — >■ oo. 
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FIG. 2. (Color online) Plot of nc,wrap (i) versus logL for model (a) on the square lattice, together with two typical fits. One fit 
(green curve) is logarithmic, nc,wrap(i) ~ 4/ ln(25I/)^'^*. It assumes that nc,wrap = nc,wrap(i ~ oo) — 0. The other fit (blue) 
assumes that ric.wrap is non-zero, and that nc,wrap(i) has power law finite size corrections: nc,wrap(i) ~ 0.031 + 0.35/1/°'^* . 
Notice that both fits involve three free parameters. The logarithmic fit is slightly better, but it is not clear how significant 
this is, since one should expect an infinite series of further correction terms in any such fit. Notice also that nc, wrap is very 
small (i.e., the average cluster mass at criticality is huge) in the second fit, in striking contrast to other (bond, site) percolation 
models on any 2-d lattice. In our opinion, this is a stronger argument against the second fit than the quality of the fit itself. 
Finally, no decent fit is possible with a pure power law nc,wrap(i) = a/L" , expected if r < 2 - although fits with two or more 
powers with different exponents would be possible. 
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FIG. 3. (Color online) Mass distributions (log-log) for model (a) on the triangular lattice with L — 4096, for three different 
values of ?i. The middle (green) curve has the longest straight piece and is thus used to define nc,T- The other two curves show 
the sub- (super-)critical behavior. For increased significance, we actually do not show P„(m) but rather m^P„(m). The fact 
that the central (green) curve has a negative slope clearly indicates r > 2. Notice that this conclusion is rather robust and 
would not require a very precise estimate of nc,T. 
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FIG. 4. (Color online) Same as Fig. S3, but for the square lattice. Superficially, the graph looks similar to Fig. S3, but the 
slope of the critical (green) curve now is positive, clearly suggesting that r < 2 unless there are very strong corrections to 
scaling. 
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FIG. 5. (Color online) Similar to Fig. 5 of the main paper, but for the fractal dimension of the largest (percolating) cluster. 
It is obtained from the positions of the peaks of the central (critical) curves in mass distribution plots such as those in Figs. 
S4 and S5. Again we see that both curves increase with L. Again the triangular lattice data are compatible with ordinary 
percolation, D = 91/48 = 1.8958 . . ., but the square lattice data are not. There, the percolating cluster has a distinctly larger 
dimension which is moving away from the ordinary percolation value. The analytical arguments given in the main paper allow 
only three possible scenarios: (i) D < 2,t > 2,v < oo, and power law corrections to scaling; (ii) D = t = 2,v = oo, and 
logarithmic corrections; and (iii) D = 2,t < 2,v < oo, and power law corrections to scaling. Scenario (i) (i.e., r > 2) seems 
definitely ruled out for the square lattice by the data. Scenarios (ii) and (iii) are both possible, but would require very large 
and slowly vanishing finite size corrections for Dcfi(L). At face value, the data for nc{L) shown in Fig. S2 and in Fig. 3 of the 
main paper would favor scenario (ii), while Tefi{L) (Fig. 5) would favor scenario (iii). Our final preference of scenario (ii) is 
based on the fact that it is the only one which can naturally accommodate logarithmic (and thus very slowly vanishing) finite 
size corrections. 



